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Halo concentrations in the standard ACDM cosmology 
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ABSTRACT 

We study the concentration of dark matter halos and its evolution in N-body sim- 
ulations of the standard ACDM cosmology. The results presented in this paper are 
based on 4 large iV-body simulations with ~ 10 billion particles each: the Millennium-I 
and II, Bolshoi, and MultiDark simulations. The MultiDark (or BigBolshoi) simula- 
tion is introduced in this paper. This suite of simulations with high mass resolution 
over a large volume allows us to compute with unprecedented accuracy the concen- 
tration over a large range of scales (about six orders of magnitude in mass), which 
constitutes the state-of-the-art of our current knowledge on this basic property of 
dark matter halos in the ACDM cosmology. We find that there is consistency among 
the different simulation data sets, despite the different codes, numerical algorithms, 
and halo/subhalo finders used in our analysis. We confirm a novel feature for halo 
concentrations at high redshifts: a flattening and upturn with increasing mass. The 
concentration c(M, z) as a function of mass and the redshift and for different cosmo- 
logical parameters shows a remarkably complex pattern. However, when expressed in 
terms of the linear rms fluctuation of the density field <j(M, z), the halo concentration 
e(cr) shows a nearly-universal simple U-shaped behaviour with a minimum at a well 
defined scale at a ~ 0.71. Yet, some small dependences with redshift and cosmol- 
ogy still remain. At the high-mass end (er < 1) the median halo kinematic profiles 
show large signatures of infall and highly radial orbits. This c-a(M, z) relation can 
be accurately parametrized and provides an analytical model for the dependence of 
concentration on halo mass. When applied to galaxy clusters, our estimates of concen- 
trations are substantially larger - by a factor up to 1.5 - than previous results from 
smaller simulations, and are in much better agreement with results of observations. 

Key words: cosmology: theory - dark matter - galaxies: halos - methods: N-body 
simulations. 



1 INTRODUCTION 

N-body cosmological simulations have been essential for 
understanding the growth of structure in the Universe, 
and in particular, they have been crucial for studying the 
properties of dark matter halos in the standard Lambda 
Cold Dark Matter (ACDM) cosmology. Simulations are 
also an invaluable tool for analyzing galaxy surveys, for 
studying the abundance evolution of clusters of galax- 
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ies, and for semi-analytical models of galaxy formation. 
In recent years, the development of numerical codes and 
access to powerful supercomputers have made it possi- 
ble to perform Grand Challenge cosmological simulations 
with high mass resolution over a large volume, which pro- 
vide the basis to attack many problems in cosmology. 
The Millennium simulation i|Springel et al.ll2005l . MS-I), its 
smaller volume and higher-r esolution Millennium-II version 
i|Bovlan-Kolchin et alfhopgl. M S-ID and the new Bolshoi 
simulation ( Klypin et al.l 201(ih . in this respect, constitute 
a remarkable achievement. 
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Millennium, Bolshoi, and the new MultiDark (or Big- 
Bolshoi) simulation, which is introduced in this work, al- 
low us to estimate for the ACDM cosmology, with unprece- 
dented statistics, the concentration of dark matter halos 
and its evolution over six orders of magnitude in mass. The 
comparison between this suite of cosmological simulations 
also provides an unique opportunity to study consistency 
on dark matter halo statistics using different codes, numer- 
ical algorithms, halo finders and cosmological parameters. 
Halo concentrations c = i? v i r /V s have been studied exten- 
sively during the last decade. Here -R v ir is the virial ra- 
dius of a halo and r s is the break radius between an in- 
ner ~ r _1 density profile and an outer r -3 profile. Typi- 
cally, the median halo concentration declin es with increas- 
ing mass and redshift (|Bullock et al.|[200ll ). and the shape 
of the mass-concent ration median relation evolves. Recently, 
iKlypin et al.l (|2010| ) found a novel feature: at high redshift, 
the concentration flattens, and then increases slightly for 
high masses. Numerous works made use of cosmological sim- 
ulations to determine the correlation and scatter of concen- 
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environment and cosmologv (e.g. Navarro. Frenk, & White 


1997; Bullock et al.l 2001; Eke. Navarro. & Steinmctz 


2001; 


Wechsler et al.l 


20021; Zhao et all 20031; Neto et al. 




2007; 


Gao et alJ I200S 


; iMaccio. Dutton. & van den Bosch 




2008; 


Zhao et al.l 2008 


; iMunoz-Cuartas et al.l 20ld; Klvpin et all 


2010TI. Moreover, the dependence of halo concentrations on 



their mer ging history ha s been s tudied also i n det a il. For 
example, Iwechsler etail (|2002h : IZhao et all (|2003l . l2009h 
determine concentrations using accurate modeling of halo 
mass accretion histories. Simulations and analytical mod- 
els have tried also to understand the formation and evo- 
lution of the central density cusp of dark matter halos - 
assuming that they grow inside out during the different 
phases in their mass accretion histories - to provide pre- 
dictions on hal o concentrations (see e.g. iReed et al.[ 2005; 
Lu et ail 120061; iRomano-Di'az et ail 120071 ; iHoffman et ail 



20071 ; ISaryador-Sole et al.ll2007l 1. The origin of halo concen- 



trations ultimately is linked to the physical process that 
yield the formation of the NFW density profiles in the cold 
dark matter scenario (see also iDalal. Lithwick. fc KuhlerJ 
120101 ). 

Here, we show that the dependence of concentration on 
halo mass and its evolution can be obtained from the rms 
fluctuation amplitude of the linear density field <r(M, z) by 
using our new analytical model based on the concentration- 
o~(M, z) relation. This model is able to reproduce all relevant 
features observed in the halo mass-concentration median 
relations at different redshifts: the decline of concentration 
with mass, its flattening and upturn. 

In this paper, we study the evolution of halo concentra- 
tions in N-body cosmological simulations of a ACDM cos- 
mology. We cover a large range of scales, going from halos 
hosting dwarf galaxies to massive galaxy clusters. This cor- 
responds to halo maximum circular velocities ranging from 
25 to 1800 km s _1 , covering six orders of magnitude in mass. 
In Section 2 we describe the simulations and halo catalogs 
used in this work. Methodology to estimate halo concen- 
trations is presented in section 3. In this paper we analyze 
only distinct halos. A halo is called distinct if its center is 
not inside the virial radius of a larger halo. In Section 4 we 
study concentrations of distinct halos and their evolution 



with redshift using data from the Millennium, Bolshoi and 
the new MultiDark simulations. A new parameterization of 
the distinct halo concentrations as a function of the rms 
of the linear density field o(M, z) is presented in Section 5. 
Main results of this work are summarized in Section 5. An 
Appendix gives discussion of various selection effects. 



2 SIMULATIONS AND HALO 
IDENTIFICATION 

Table 1 summarizes the basic numerical and cosmological 
parameters of the four simulations used in this work. The 
Millennium simulation dSpringel et ail 20051) and the more 
recen t Millennium-II simulation ( Bovlan-Kolchin et ail 
2009) adopted the same WMAP first year cosmological pa- 
rameters, used the same number of particles to resolve the 
density field, and share the same output data structure. MS- 
II was done in a cubic simulation box one-fifth the linear size 
of the original MS-I with 5 times better force resolution and 
125 times better mass resolution. MS-I, II were run with the 
TREE-PM codes GADGET-2 and GA DGET-3 respectively. 

The Bolshoi (|Klypin et ail 120101 ) and the MultiDark 
simulations use the latest WMAP5 and WMAP7 cosmo- 
logical parameters, which are also consistent with other re- 
cent o bservational constraints - see Figure 1 of lKlypin et ail 
(|2010h . The Adaptive-Refinement- Tree (ART) code was 
used for these two simulations. ART is an Adaptive-Mesh- 
Refinement (AMR) type code. A d e tailed desc ription of 
the co de is given in iKravtsov et ail |l997) and iKravtsovl 
(1999). The code w as parallelized using MPI li braries and 
OpenMP directives l|Gottloeber fc KlvpinlbOOST ). Details of 
the time-stepping algorithm and c omparison with GAD GET 
and PKDGRAV codes are given in lKlypin et~ail (|2009h . The 
ART code increases the force resolution by splitting individ- 
ual cubic cells into 2x2x2 cells with each new cell having 
one half the size of its parent. This is done for every cell 
if the density of the cell exceeds some specified threshold. 
The value of the threshold varies with the level of refinement 
and with the redshift and is typically 2-5 particles per cell. 
180 snapshots were save d for analysis. Detai ls of the Bolshoi 
simulations are given in IKlypin et all i|2010r i. 

Here we give some details of the MultiDark run. Initial 
conditions were set at the redshift z- m it = 65 using the same 
power spectrum as for Bolshoi. The force resolution in the 
ART code varies with time. The comoving resolution (size of 
the smallest cell) for the MultiDark was equal to 244 /i _1 kpc 
at z > 10. When the fluctuations started to collapse, the 
resolution became smaller and at z = it was 7.6/i _1 kpc. 
The ART code is designed in such a way that the proper 
(physical) resolution is nearly preserved over time. For Mul- 
tiDark the proper resolution was ~ 7/i _1 kpc for z — — 8 
as compared with l/i _1 kpc for z — — 20 for Bolshoi. 
The total volume resolved at the highest resolution at z = 
was 360(/i _1 Mpc) 3 and 2.8 x 10 4 (/i _1 Mpc) 3 for twice worse 
resolution. 50 snapshots were saved for analysis. 

The main difference in cosmological parameters be- 
tween the simulations is that the Millennium simulations 
adopted a substantially larger amplitude of perturbations 
us that is nearly 4cr away from r ecent constraints (e.g., 
IKlypin et ail 120101 : IGuq et alj|201ll ). The difference is even 
larger on galaxy scales since the Millennium simulations 
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used a large r tilt n 3 of the powe r spectra (see Table 1, and 
Figure 2 in iKlypin et all |2010l i for the comparison of the 
linear power spectra of Bolshoi and Millennium simulations). 

In this paper we analyze only distinct halos. A halo is 
called distinct if its center is not inside the virial radius of a 
larger halo. One of the most important characteristics of a 
halo is its maximum circular velocity: 



Vmax = max 



GM(< r) 



(1) 



There are several advantages of using V max to character- 
ize a halo as opposed to the virial mass. First, V ma x does 
not have the ambiguity related with the definition of "virial 
mass." Virial mass and radius vary depending on overdensity 
threshold used. For the often-employed overdensity 200 and 
"virial" overdensity thresholds, the differences in definitions 
result in changes in the halo radius from one definition to an- 
other and, thus, in concentration, by a factor of 1.2-1.3, with 
the exact value being dependent on the halo concentration. 
Second and more important, the maximum circular velocity 
Knax is a better quantity to characterize halos when we re- 
late them to the galaxies inside these halos. For galaxy-size 
halos the maximum circular velocity is defined at a radius of 
~ 40 kpc: closer to sizes of luminous parts of galaxies than 
the much larger vi rial radius, which fo r the Milky- Way halo 
is ~ 250 kpc fe.g.. lKlypin et alj|2003 ). 

In both Millennium simulations dark matter halos 
were found using the friends-of-friends algorithm (FOF, 
iDavis et all 1 19851). with a linking len gth of b = 0.2. The 
SUBFIND code l|Springel et alj 120011 ) searched for bound 
substructures within every FOF group. In most cases, the 
main component of a FOF group is typically a dominant 
spherical "subhalo" whose particles make up most of the 
bound component of the FOF group. This dominant or cen- 
tral subhalo can be identified as a distinct halo with a maxi- 
mum circular velocity Vmax- For our analysis we consider all 
distinct halos (i.e. dominant subhalos of FOF groups) with 
maximum circular velocities greater than Vm ax > 25 km s _1 
and Vmax > 70 km s _1 for MS-II and MS-I respectively. 
Halo velocity functions a r e complete above these velo cities 
(see ISpringel et ahl 120051 ; iBovlan-Kolchin et al l I2009D . We 
have downloadecQ MS-I and MS-II halo catalogs for several 
epochs up to z — 10, as close as possible to those redshifts 
available for Bolshoi and MultiDark. 

Dark matter halos are identified in Bolshoi and Multi- 
Dark with a paralle l version of the Bound-D ensity-Maxima 
(BDM) algorithm jKlypin fc Holtzmadfl997l ). The BDM is 
a Spherical Overdensity (SO) code. It finds all density max- 
ima in the distribution of particles using a top-hat filter 
with 20 particles. For each maximum the code estimates the 
radius within which the overdensity has a specified value. 
Among all overlapping density maxima the code finds one 
that has the deepest gravitational potential. The position of 
this maximum is the center of a distinct halc0. 

The halo radius can be defined as the radius of a sphere 
within which the average density is A times larger than the 



1 FOF and SUBFIND halo catalogs from MS-I a nd MS-II are 
available at http://www.g-vo.org/MyMillennium3 . 

2 Halo catalogs for these simulations are available at 
http://www.multidark.org 
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Figure 1. Mass functions of distinct halos in Bolshoi (open cir- 
cles) and MultiDark (filled circles) simu l ation s at z = 0. The 
analytical approximation of lTinker et al.l (2008) is shown by the 
curve. In the overlapping region both simulations produce nearly 
identical results. 



critical density of the Universe 
4ty 



M 



Apcr (Z)-R 



(2) 



where p C r(z) is the critical density of the Universe. In this 
paper we use the overdensity A = 200 threshold. Corre- 
sponding values for mass and radius are M = M200 and 

R = i?200- 

Figure[T]shows differential mass functions of distinct ha- 
los in Bolshoi and MultiDark simulations. The two simula- 
tions demonstrate remarkably good consistency in the over- 
lapping region of masses (10 12 - 1O 14 )/i _1 M . To describe 
the mass function dn/dM we use the following equations: 



dn 

1m 



TI7 = /(-) 



= A 



p m d log a 
M dM 
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ex P — 



a 2 (M, 



1 

2^2 



P(k,z)W (k,M)k 2 dk, 



(3) 
(4) 
(5) 



where p m is the cosmological matter density, and P(k, z) 
is the power spectrum of fluctuations. Here a(M,z) is the 
linear rms fluctuation of density field on the scale M as 
estimated with the top-hat filter, whose Fourier spectrum is 
W(k,M). Parameters o f the a pproximation are nearly the 
same as m iTinker et all (^OOSt ) for the same overdensity A: 
A = 0.213, a = 1.80, b = 1.85, and c = 1.57. 

Analysis of halo concentrations is done for different pop- 
ulations of halos. We limit our analysis to halos containing 
more than 500 particles. Most of the time we analyze all 
such halos regardless of their degree of virialization or in- 
teraction. Occasionally we select only "relaxed" halos. A 
number of parameters is involved in this selection. We de- 
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Table 1. Basic parameters of the cosmological simulations. L box is the side length of the simulation box, N p is the number of 
simulation particles, e is the force resolution in comoving coordinates, M p refers to the mass of each simulation particle, and the 
parameters f2 m , £7 a , Qf,, n s (the spectral index of the primordial power spectrum), h (the Hubble constant at present in units 
of 100 km/sMpc' 1 ) and <x 8 (the rms amplitude of linear mass fluctuations in spheres of 8 h 1 Mpc comoving radius at redshift 
z = 0) are the ACDM cosmological parameters assumed in each simulation. 



Name 


^box 

(/i^Mpc) 


N p 


e 

(fc-ikpc) 


M p 




n A 


n b 


n s 


h 


CT8 


reference 


Millennium 


500 


2160 3 


5 


8.61 10 s 


0.25 


0.75 


0.0450 


1.00 


0.73 


0.90 


Springel et al. 2005 


Millcnnium-II 


100 


2160 3 
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6.89 10 6 
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0.75 


0.0450 


1.00 


0.73 


0.90 


Boylan-Kolchin 
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2048 3 


1 
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0.82 


Klypin et al. 2010 
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this work 
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Figure 3. Examples of halo density profiles (full curves) and their fits (dashed curves) using concentrations obtained with the ratios 
of the maximum circular velocity Vinax to the halo velocity V200- Each panel shows two full curves: the density profile of all particles 
(upper curve) and only bound particles (lower curve). Vertical dotted lines show the outer radius of bound particles. Typically the 
highly concentrated and relaxed halos are well fit by the approximations (two bottom panels). Low concentrated halos and halos with 
large perturbations (two top panels) show large deviations between the data and approximations. Fitting those unrelaxed profiles with 
standard minimization technique does not improve the accuracy of the fits as indicated by dotted curve in the top-right panel. 
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Figure 2. Relation between halo concentration c and the ratio 
of the maximum circular velocity Vmax to the circular velocity 
at the virial radius V200 f° r the NFW profile (full curve) and for 
three Einasto profiles with parameters a = 0.15,0.20,0.25 (from 
top to bottom). 



fine an offset parameter x g as the distance between the 
center of a halo (density maximum with the deepest grav- 
itational potential) and the center of mass of the halo in 
units of the virial radius. For each halo we also use the spin 
parameter A and the virial ratio 2K/\U\ — 1, where K and 
U are the total kinetic and potential energies. A halo is 
called "relaxed" if it satisfies the followin g three conditions : 
i„ff < 0.1, 2K/\U\ - 1 < 0.5, and A < l.lNeto et all (|2007T I 
and iMaccio. Dutton. fc van den Boschl (|2008l ) give detailed 
discussion of these conditions. 



3 FINDING HALO CONCENTRATIONS 

There are different ways to estimate halo concentration. Tra- 
ditionally it is done first by fitting a NFW profile 

Ap a _ r 



Pnfw 



(6) 



x(l + x) 2 ' r 3 

to the spherically averaged density profile of a halo, where 
p s and r 3 are the characteristic density and radius. Then the 
halo concentration is found as the ratio of the virial radius 
to the characteristic radius r 3 : 

c = (7 ) 
r a 

However, the fitting can be difficult. For example, special 
attention should be paid to the central region: if fitting 
starts too close to the center where the resolution is not 
sufficient, the fitting will produce too low concentration be- 
cause the density in the center is underestimated. In addi- 
tion to the problems with the center and the necessity to 
deal with binned data, direct fitting has another feature: it 
heavily relies on an assumption of a particular shape of the 
halo profile, which is typically taken as the NFW profile. 



Halo profiles on average show systematic deviations from 
the NFW approxima tion and are better approximated with 
the Einasto profile (|Navarro et alj 120041 : IGao et all 120081 : 
iNavarro et al.|[2010l ). One may naively expect that fitting a 
more accurate approximation would produce better results 
on concentrations. This is not the case. There is a reason 
why halo concentrations are not estimated with the Einasto 
profile: results show large fluctuations in the concentration, 
which are due to a smaller curvature of density profile p(r) in 
the crucial range of radii around the radius with the logarith- 
mic slope -2. These weaknesses of the direct fitting provide 
a motivation to look for alternatives. 

The idea is to start with the NFW profile. It is uniquely 
defined by the characteristic density p a and the character- 
istic radius r s . Instead of these two parameters one can use 
any two independent quantities to define the profile. For 
example, one can choose the virial radius and the radius at 
whic h the circular veloc ity reaches half of its maximum value 
as m lAlam et all (|2002T l. Or one can use the radius contain- 
ing 1 /5 of virial mass and the virial mass |Avila-Reese et all 
1999). Recently iKlypin et al.l (|201Ch proposed to use a rela- 
tion between the maximum circular velocity Knax, the virial 
mass and the virial radius to estimate the halo concentra- 
tion. We follow this idea and use a more simple version with 
parameters Kn ax and the circular velocity at the virial ra- 
dius: 



V200 



/ GM200 

V ^200 



1/2 



(8) 



In this case the ratio Knax/Vioo of the maximum circular 
velocity Vmax to the virial velocity V200 is directly related 
with the halo concentration and can easily be interpreted 
because larger ratios imply larger concentration. The value 
of the halo concentration can be estimated by assuming a 
shape of the density profile. However, this is not required 
because the velocity ratio is a measure of the halo concen- 
tration regardless of what profiles halos have. For the case of 
the NFW halo density profile, the Vmax/V-soo velocity ratio 
is given by the following relation: 



V200 



0.216c 

1W 



1/2 



where /(c) is 



/(c) = ln(l + c) 



(9) 



(10) 



Having the Vmax/Vioo ratio for each halo, we find the halo 
concentration c by solving numerically eqs. - I10p . For 
convenience, Figure [2] shows the dependence of concentra- 
tion c on Vmax/V^oo in a graphical form . We also present re- 
sults for three typical Einasto profil es (|Navarro et al ] 12004 
IGao et aLlliooi : INavarro et afll2010h : 



PEin = Ps exp 



(x a - 1) 



r 



(11) 



In practice we bin halos by their maximum circular ve- 
locity Vmax- At a given redshift and for each circular velocity 
bin we obtain the median virial velocity V^oo- Then we com- 
pute c by solving eqs. (P-fTO)). 

We made different tests of the accuracy of our algo- 
rithm. We find quite a good agreement of our es timates of c 
with the results of direct fitting of halo profiles bv lNeto et al.l 




Figure 4. Comparison of halo concentrations cat obtained by 
fitting the NFW profiles to halo densities with concentrations 
Cyratio found using the ratios of the maximum circular velocity 
Vmax to the halo velocity V200- The most massive ~ 10000 halos 
in the Bolshoi simulation are used for the plot. The top (bottom) 
panel shows halos at z = (z = 3). All halos are used regardless 
of their relaxation status. Full lines present the relation cv ra tio = 
c NFWflt- The dashed lines show a 15% offset between the two 
concentrations. Most of the median concentrations fall between 
those two lines. Error bars present 10 and 90 percent spread of 
the distribution of individual concentrations. 



l|2007l ) with the deviations between the two methods being 
less than 5 percent for the whole range of masses in the MS 
simulation. Figure [3] shows six examples of fitting density 
profiles of large halos with M 200 « (1.5 - 3) x 10 14 /i _1 M© 
in the Bolshoi simulation. Each of the halos is resolved with 
more than a million particles and the random noise associ- 
ated with finite number of particles is practically negligible. 
We use 30 radial bins equally spaced in logarithm of radius 
with A log(r/J?20o) = 0.05 up to R — 2i?200- Plots show typ- 
ical examples of high and low concentration halos. Just as 
expected, for halos with smooth density profiles (the bottom 
row of panels) the accuracy of fits obtained with Vmax/V^oo 
concentrations is very good. Halos with large substructures 
and big radial fluctuations (typically halos with low con- 
centration; the top row of panels) result in relatively poor 
fits. However, in these cases the direct fitting does not pro- 
duce much better approximations as demonstrated by the 
top right panel. 

Figure [4] provides a more extensive comparison of con- 
centrations obtained with the direct fitting and with the 
Knax/V^oo ratio method. For this plot we select all ~ 7000 
halos in the Bolshoi simulation with Vm ax > 350 km s _1 at 
z = and ~ 10000 halos with Vmax > 250 km s" 1 at z = 3. 
The halos with concentration larger than c w 5 demonstrate 
a tight relation between the two estimates of concentra- 
tion: more than 90 percent of halos fall between two lines: 
cvratio = est and cvratio = 1.15cflf At smaller concentrations 



Figure 5. The ratio Vmax/Vaoo of the maximum circular velocity 
to the virial velocity as a function of mass M200 for distinct halos 
at different redshifts for MS-I (filled symbols) and MS-II (open 
symbols) simulations. Error bars are statistical uncertainties. The 
MS-I and MS-II simulations agree quite well at z = 0. At higher 
redshifts there are noticeable differences between MS-I and MS-II. 

there is a small systematic (5-15)% offset between the two 
estimates. 

These results demonstrate that both the direct fitting 
of density profiles and the method based on the V max /V2oo 
ratio give similar concentrations. There are the small sys- 
tematic deviations, which for relaxed halos or for halos with 
concentrations c > 5 are less than a few percent. For halos 
with small concentrations the deviations are larger, but are 
still less than ~ 15%. We ignore these differences and use 
concentrations based on the I^nax/Vioo ratios. 



4 EVOLUTION OF THE DISTINCT HALO 
CONCENTRATIONS WITH REDSHIFT: 
RESULTS OF SIMULATIONS 

Our analysis of halo concentrations in different simulations 
starts with the simpler statistics of Vmax/V2oo- It is directly 
related with the concentration c, but V max /V2oo does not 
require an assumption of a particular density profile. These 
raw data are presented in Figures [5]|6] for all simulations 
discussed in section 2. Here we consider only halos with more 
than 500 particles for our analysis. Error bars correspond 
to the statistical uncertainty in the determination of the 
Vmax/V2oo median velocity ratio. Halos are binned by Vmax, 
and the average M200 for each velocity bin is the horizontal 
axis. 

At 2 = 0, the agreement between both Millennium sim- 
ulations is very good in the range of masses where they 
overlap. This is expected if the two simulations faithfully 
present results of the cosmological model and are not af- 
fected by finite box sizes or any other numerical issues. Yet, 
at higher redshifts there are notable differences between MS- 
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Figure 6. The same as Figure [5] but for Bolshoi (open sym- 
bols) and MultiDark (filled symbols) simulations. Both simula- 
tions show remarkable agreement at all masses and redshifts. 



I and MS-II (see Figure [5}. At a fixed halo mass M200, data 
from MS-I yield systematically larger V ma x/V2oo velocity ra- 
tios, increasing with redshift up to a 3 — 4% difference. For 
example, at z = 3 we find ~ 3% difference for 10 12 M©. 
This will translate to a difference of about 20 — 40% in terms 
of halo concentration according to eq. ©. 

Regardless of this discrepancy between the two Millen- 
nium simulations, we do clearly see an upturn of the velocity 
ratios, and, hence, halo concentrations with increasing halo 
mass at higher redshifts. The same is observed in the Bolshoi 
and MultiDark simulations (see Figure [6]). We conclude that 
concentrations of distinct halos do increase w ith halo mass 
at hig her redshifts as previously reported by iKlypin et al.l 
(|2010f ) . In spite of the qualitative agreement, there are differ- 
ences between Bolshoi/MultiDark and the Millennium-I/II 
simulations. Those are expected and are due mainly to the 
different cosmological parameters used in the simulations. 

The reason for the discrepancy between Millennium 
simulations found in the Knax/V^oo velocity ratios at high 
redshifts is not clear. It might be due to differences in mass 
and force resolution. If that were true, one may expect that 
a lower resolution MS-I simulation should have a lower con- 
centration. However, this is not the case: MS-I has larger 
concentration for overlapping masses. An important aspect 
that could cause this discrepancy is the fact that substruc- 
ture within the dominant (distinct) halo is not taken into 
account when Kn ax is computed. Certainly this could make 
differences up to the level of few percent as observed. The 
MultiDark simulation is helping us to clarify this issue be- 
cause MultiDark and Bolshoi have very different resolutions. 
In this case, we do see in Figure [5] a good agreement with 
the Bolshoi simulation at all redshifts in the mass interval 
where both simulations overlap. The BDM halo finder does 
include substructures within distinct halos to compute V~ ma x. 

So far we were following the traditional path of studying 
concentrations by expressing results for c as functions of halo 




Figure 7. The ratio F ma x/V2oo as a function of log a~ 1 . In this 
plot halos mass increases from left to right. The large variations 
in shapes and amplitudes of F m ax/V200 — AT seen in Figures 1 5161 
are replaced with much tighter relations Fmax/V^oo ~~ l°g c — 1 . 
All simulations show the same pattern: Fmax/V^oo has a U-shape 
with an upturn at large masses. 



mass at different redshifts c = c(M, z) ( e.g., Bullock et al.l 
200 it lEke. Navarro, fc Steinmetzl I2001I; IZhao et alJ 120031 : 



Maccio. Dutton. fc van den Boschl 12008) . As many other 



groups find, results for c(M, z) appear to be not simple: 
the concentration evolves in a complicated way. In addition 
to previously found effects, at high redshifts the amplitude 
starts to increase, which makes the whole c{M, 2) function 
very complex. 

One of reasons why c(M, z) may look too complicated 
is because "wrong" physical quantities M and z are used. 
The situation may be compared to the studies of the halo 
mass function. Once a more physi c al variable o~(M, z) was 



used (see e.g. Ijenkins et all 120011; ISheth fc Tormenl 120021 : 
IWarren et all |2006| ; iTinker et al.ll2008l ). the results became 



more transparent. In turn, this simplification allowed the 
development of accurate approximations such as that pre- 
sented by eqs. (|3I5|I . We will follow the same direction when 
we study the halo concentration. In particular, as the main 
variable we choose log a~ 1 (M, z) - a quantity that is widely 
used in the mass function analysis. Using this variable helps 
to remove most of the evolution of the concentration with 
redshift, helps to understand the differences due to the cos- 
mological parameters and the matter power spectrum, and 
allows a more clear comparison between different simula- 
tions. 

Figure [7] shows the V m ax/V2oo velocity ratio as a func- 
tion of log <7 _1 of distinct halos at different redshifts for both 
Millennium simulations (left panels) and Bolshoi/MultiDark 
(right panels). Large volume simulations MS-I & MultiDark 
provide good statistics at the high-mass end (log a -1 > 0) 
and the smaller boxes, with much higher numerical reso- 
lution (MS-II and Bolshoi), help to reach small values of 
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Figure 8. Conccntration-o-(M) median relation for distinct halos 
in the Bolshoi (open symbols) and MultiDark (filled symbols) 
simulations for different redshifts. Solid lines are predictions from 
our analytical model given by eqs. (I14t I17t . The data at different 
redshifts look very similar, but actually have a slight evolution: 
the minimum of c becomes smaller and it is reached at smaller 
log cr~ 1 

log cr -1 < (small masses). Figure [7] also shows that when 
expressed in terms of the variable log a~ , the concentra- 
tion becomes a well-behaved function with an U shape: at 
low log cr -1 it decreases with log cr -1 . Then it reaches a min- 
imum at a well defined scale log cr" 1 ~ 0.15 (cr ~ 0.707), and 
at larger log cr" 1 it begins to increase. 

In spite of the fact that the rescaling from M to log cr -1 
significantly reduces deviations in V ma x/V2oo (and, thus, in 
concentration), it does not completely remove all the devi- 
ations. There is still some residual dependence on cosmol- 
ogy: there are differences between the Millennium and the 
Bolshoi/MultiDark simulations. In addition, there is a weak 
evolution with redshift, which is more pronounced in high 
resolution runs Bolshoi and MS-II 

As discussed above, at high redshifts there is a discrep- 
ancy between MS-I and MS-II simulations resulting in a 
20 — 40% difference in concentrations for a given halo mass 
M between the two simulations. At the same time, the over- 
lap between Bolshoi and MultiDark is good at all redshifts. 
This is the reason we use Bolshoi and MultiDark to study the 
weak evolution of the c — log <r _1 relation with the redshift. 
The time evolution of the concentration can be described as 
a decrease of the minimum value of the U-shaped concen- 
tration curves with increasing redshift (from 5.1 at z = 
down to 3.7 at z = 4), and a slight shift of the position of 
the minimum to smaller values of log a -1 , i.e. from 0.2 at 
z = to 0.05 at z = 4. Figure [8] shows details of this weak 
evolution. 

The nature of the upturn in the halo concentrations 
is an interesting problem on its own. We do not make de- 
tailed analysis of halos in the upturn, but we test some sim- 
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Figure 9. Average radial velocity (V r ) (top panels) and velocity 
anisotropy /3 (bottom panels) profiles for distinct halos selected 
from the Bolshoi simulation in different bins of log cr - 1 and red- 
shifts indicated in the plots. Halos with large log (most mas- 
sive halos) clearly show preferentially radial orbits and overall 
infall pattern at large distances from the halo center. Less mas- 
sive halos (log <r — 1 < 0) on average do not have infall velocities 
and their orbits in the peripheral regions tend to be isotropic 
while retaining slightly radial velocities in the central regions. 

pie hypotheses. One may think that the upturn is due to 
non-equilibrium effects. Indeed, the halos in the upturn are 
the largest halos at any given moment, and those halos are 
known to grow very fast. Results presented in the Appendix 
clearly show that out-of-equilibrium effects do not provide 
an explanation for the upturn: selecting relaxed halos only 
increases the magnitude of the upturn. However, the fast 
growth of the halos in the upturn does change the structure 
of the most massive halos (even those that are "relaxed"). 
Figure [9] shows average radial velocity profiles (V r ) and 
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velocity anisotropy /3 profiles for distinct halos in the Bol- 
shoi simulation for different bins of log a -1 . These statistics 
can be used to indicate the dynamical state of halos that lie 
at different locations in the V ma x./V200~o'{M) median rela- 
tions displayed in Figure [7] Average halos in the high-mass 
end (i.e. log tr -1 > 0) clearly show signatures of infall given 
by large negative radial velocities. This is more clearly seen 
when expressed in terms of (V r ), but is also present in the 
velocity anisotropy f3 — 1 — of /2<t, 2 , which tends to be larger 
at large values of log cr _1 . These results indicate that orbits 
are preferentially radial for halos in the upturn part of the 
c — M relation. The situation changes for halos with smaller 
log o . They do not show an infall pattern and their ve- 
locity anisotropy is significantly smaller. Density profiles are 
less sensitive to these effects. The halos are "normal": on av- 
erage they are well described by NFW for all log a" 1 bins. 
The only indication that something is not normal is the large 
halo concentration. 
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5 THE C cr(Af) RELATION: AN ANALYTICAL 
MODEL FOR THE HALO 
MASS-CONCENTRATION RELATION 

It is then desirable to provide for the ACDM cosmology an 
approximation that describes the dependence of halo con- 
centration c on halo mass and redshift. However, approxi- 
mations which use mass and redshift as the main variables 
are prone to a severe problem: lack of scalability. An ap- 
proximation found for one set of cosmological parameters 
and for one particular redshift is not applicable to another 
cosmology or redshift. Just as with the mass function, we 
find a dramatic improvement of the accuracy of the approx- 
imations for the halo concentration once we use <j(M,z). 
Following the same line of ideas, we use a new "time vari- 
able" x defined below. The motivation comes from the linear 
growth rate D of fluctuations in the ACDM cosmology. If it 
is normalized to be unity at z — 0, then 



D(a) 



5 /n m ,o\ i/3 vttx^ 



2 V J1a, o , 
/ fU,o x 1/3 



c 3 ' 2 dx 



r< -' J [1 + x 3 ] 3 / 2 ' 

(i+*r\ 



(12) 
(13) 



where f2 m ,o and Qa,o are the matter and cosmological con- 
stant density contributions at z = 0. Note that when D is 
written in this form, there is no explicit dependence on the 
redshift z: the time dependence goes through the variable x. 

If the concentration c depended only on a, we would 
had a universal function c(cr). Similar to the mass function, 
we find that concentration is not universal and has some de- 
pendance on the redshift. We attribute this residual redshift 
dependence of c(a, a) mostly to the change in the growth rate 
of fluctuations D(a) related to the change of matter density 
parameter f2 m (a). Another possible factor is the change with 
the redshift of the slope of the power spectrum P(k) at a 
given a. We find that corrections due to the growth rate 
factor produce sufficiently accurate fits and we do not study 
here the effect due to change in the P(k) slope. 

Because all our approximations use only a and x, we 
expect that they are applicable also for different redshifts 
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Figure 10. Dependence of the halo minimum concentration c m i n 
(bottom panel) and the value of cr"^ (top panel) on x for distinct 
halos in the Bolshoi and MultiDark simulations. Solid lines are 
the approximations given in eqs. H19H20> . 



and different cosmological parameters. Fitting of all avail- 
able data from Bolshoi and MultiDark simulations gives the 
following approximation for halo concentration: 



c(M, z) 

! 

(7 

C(a') 
where 



B (x)C(a'), 
Bx(x) o(M,x), 



= A 



+ 1 



exp 



(14) 
(15) 

(16) 



A = 2.881, b = 1.257, c = 1.022, d = 0.060 (17) 

The parameterization form of the function w as motivated 
by the Sheth-Tormen jSheth fc Tormenl l2002h approxima- 
tion for the halo mass function. Here, the functions Bq(x) 
and Bi(x) are defined in such a way that they are equal to 
unity at z = for the WMAP-5 parameters of the Bolshoi 
and MultiDark simulations. Thus, the function C(cr') is the 
concentration at z — for this cosmological model. We find 
the following approximations for Bq(x) and Bi(x): 



B (x) 



ain (^) 



c min (1.393)' 



Bi(x) 



"min (1-393)' 



(18) 



where c m i n and define the minimum of the halo con- 

centrations and the value of a at the minimum: 



Cmin(£) = C + (Cl - C ) 

ownO) = fo" 1 + or 1 - &o l ) 

where 

c = 3.681, ci = 5.033, a = 6.948, x = 0.424, 



— arctan [a(x — xq)] + — 
tt 2 

— arctan \8(x — xi)] + ^ 

7T 2 



(19) 
(20) 



(21) 
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Figure 11. Dependence of halo concentration c on log a~ 1 after 
rescaling all the results of Bolshoi and MultiDark simulations to 
z = 0. The plot shows a tight intrinsic correlation of C on a' . 

and 

ao 1 = 1.047, o-jT 1 = 1.646, /3 = 7.386, an = 0.526. (22) 

Accurate approximations for the rms density fluctua- 
tion a(M,a) for the cosmological par ameters of the Bol - 
shoi/MultiDark simulations are given in lKlypin et alj (|2010T l 
and for convenience are reproduced here: 

16 9 v 0A1 

cr(Af , a) = D(a) l+ J + M , (23) 



[W 12 h- 1 M Q _ 

Figure [TO] shows the evolution of c m i n and cr~ in with 
"time" x, and presents the approximations given in eqs.()19t- 
I20p . The evolution is clearly related with the transition from 
the matter dominated period (Sl m (a) w 1, x < 0.3) to the 
A-dominated one with x > 0.7. Approximations for the 
halo concentration are presented in Figure [5] for some red- 
shifts. The parameters A, b, c, d of the C(a') relation are de- 
termined from the best fit to the concentration-cr(Af) Bol- 
shoi/MultiDark data at all redshifts. 

Here is a step-by-step description how to estimate halo 
concentration: 

• For given mass M and a = 1/(1 + z) find x, D(a), and 
a(M,a) using eqs. (fT3l [T2l [T6l or 

• Use eq. (|18p to find parameters Bo and B\. 

• Use eqs. (|15I16[1 to find a' and C 

• Use eq. (|14[) to find halo concentration c(M, z). 

We present the final results and approximations in two 
different forms. Functions Bo and B\ can be used to find 
values of C and a', which is effectively the same as rescal- 
ing concentrations c(cr, x) measured in simulations to the 
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Figure 12. Halo mass— concentration relation of distinct halos at 
different redshifts in the Bolshoi (open symbols) and MultiDark 
(filled symbols) simulations is compared with analytical approxi- 
mation eqs. (ll4ll6l (curves)). The errors of the approximation are 
less than a few percent. 

same redshift z — 0. Figure [TT] shows results of simulations 
rescaled in this way. The U-shape of C(cr') is clearly seen. 
The C(a') function to some degree plays the same role for 
concentrations as the function f(a) for the mass function 
in eqs. (|3l4fl . It tells us that there is little evolution in the 
dependence of concentration with mass once intrinsic sear- 
ings (e.g., x instead of expansion parameter) are taken into 
account. 

Another way of showing the approximations is simply 
plot eqs. (|14I16[I for different redshifts and compare the re- 
sults with the median concentration - mass relation in our 
simulations. This comparison is presented in Figure [T3I It 
shows that the errors of the approximation are just a few 
percent for the whole span of masses and redshifts. 



6 SUMMARY AND CONCLUSIONS 

We study the halo concentrations in the ACDM cosmology, 
from the present up to redshift ten, over a large range of 
scales going from halos similar to those hosting dwarf galax- 
ies to massive galaxy clusters, i.e. halo maximum circular 
velocities ranging from 25 to 1800 km s _1 (about six orders 
of magnitude in mass), using cosmological simulations with 
high mass resolution over a large volume. The results pre- 
sented in this paper are based on the Bolshoi, MultiDark, 
and Millennium-I and II simulations. There is a good con- 
sistency among the different simulation data sets despite 
the different codes, numerical algorithms, and halo/subhalo 
finders used in our analysis. 

The approximations given here for the evolution of the 
halo concentration constitute the state-of-the-art of our cur- 
rent knowledge of this basic property of dark matter halos 
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Figure 13. Comparison of observed cluster concentrations (data 
points with error bars) with the prediction of our model for me- 
dian halo concentration of cluster-size halos (full curve). Dot- 
ted lines show 10% and 90% percentiles. Open circles show re- 
sults for X-ray luminous galaxy cluste rs observed with X MM- 
Newton in the redshift range 0.1-0.3 (Ettori et alJl201Ch . The 
pe ntagon presents g a laxy k inematic estimate for relaxed clusters 
by^ 
by 



Woitak fe Lokad <201Ch. The dashed cu rve shows prediction 



Maccio, Dutton, fc van den Boschl J2008h . which significantly 



underestimates the concentrations of clusters. 



found in N-body ACDM simulations. Naturally, new simu- 
lations that improve mass resolution in even larger volumes 
are needed to face the challenges imposed by current and fu- 
ture observational programs. Our analysis can also be useful 
for comparison with analytical works that aim to understand 
the statistics and structural properties of dark matter halos 
in the standard ACDM cosmology. 

It is interesting to compare these results with other 
simu lations and models in the literature. IZhao et alj l|2003l. 
I2009D were the first to find that the concentration flattens 
at large masses and at high redshifts. Actually, figure 15 
m IZhao et al l (|2009h also shows an upturn in concentra- 
tion at z — 4. However, the authors do not even mention it 
in the text. In addition, their model prediction of the halo 
mass-concentration relation, based on the halo mass accre- 
tion histories, failed to reproduce the upturn behavior of the 
concentration with increasing mass. 

For M vir = 10 12 /t~ 1 M(7 ) in the M S-II and Aquarius sim- 
ulations iBovlan-Kolchin et alj (|2009l ) give concentration of 
c v i r = 12.9, which is 1.3 times larger than what we get from 
Bolshoi. Most of the differences are likely due to the larger 
amplitude of cosmological fluctuations in MS simulations be- 
cause of the combination of a larger ag and a steepe r spec - 
trum of fluctuations. On larger masses iNeto et all (|2007T I 
give the following approximation for all halos for MS cos- 
mological parameters: c 2 oo = 7.75(M2oo/10 12 /i" 1 M Q ) _(1 - 11 . 
Thus, the MS-I has a small (~10%) difference in C200 as 
compared with our results for M v j r = 10 14 — lO 15 /i _1 M0. 

If we use the same selection conditions (all ha- 



los selecte d by mass) and use the same cosmo logical 
model as in lMaccio. Dutton. fc van den Boschl l|200ct ). then 
the concentr ations of Milky- Way-size halos (M200 = 
10 12 h~ 1 M Q in lMaccio. Dutton. fc van den Boschl i|200Sl ) are 
10% lower than what we find - a reasonable agree- 
m ent. However, our results are in contradic tion with those 
of iMaccio. Dutton. fc van den Boschl {2008) when we con- 
sider clusters of galaxies. For example, for M200 = 5 x 
10 14 /i~ 1 M(7 ) we find C200 = 4.6 while a pproximations in 
M accio. Dutton. fc van den Boschl (|2008l ) give a substan- 
tially lower value of C200 = 3.1 - a 50% smaller value. 

In Figure [T3] we compare predictions of our model 
eqs. (|14fll6|) for the median mass-concentration relation for 
galaxy clusters (solid line) with observational estimates 
obtained from sta te-of-the-art esti mates based on X-ray 
i|Ettori et all l201Ch and kinematic l|Woitak fc Lokasl l20ldh 
data. The model has been computed at redshift z — 0.15, 
the median of the redshift distribution of the X-ray sam- 
ple. Open circles are C-M200 measurements for the sample 
of 23 X-ray luminous galaxy clusters observed with XMM- 
N ewton in the redshi ft range 0.1-0.3 selected from Table 2 
111 lEttori et all (|2010l ). We adopted for each cluster an av- 
erage concentration obtained from both concentration mea- 
surements using two different techniques applied by the au- 
thors to recover the gas and dark matter profiles. We re- 
jected those clusters where their concentration estimates 
differs more than 30% between both mass reconstruction 
techniques. The mass-concentration mean estimate from the 
kinematic analysis of the combined sample of 41 nearby clus- 
ters is also shown as a pentag on symbol (IWoitak fc Lokasl 
|2010| ). The original estimate of IWoitak fc Lokasl l|201Ch has 
been corrected to our overdensity definition. Our model pre- 
diction is in good agreement with the mass-concentration 
observational measurements. Yet, the median halo mass- 
concentration relation from Maccio et al. 2008 (dashed line) 
produces by far lower estimates of the concentration for a 
given value of the halo mass as compared to the observa- 
tional data. 

Our main results can be summarized as follows: 

• Our study of the evolution of the Vmax/V^oo veloc- 
ity ratio (which is a measure of concentration) as a func- 
tion of halo mass M200 has confirmed that the halo mass- 
concentration relation shows a novel feature at high red- 
shifts: a flattening and an upturn at the high-mass end. 

• When expressed in terms of the variable log cr _1 , the 
halo concentration is a well-behaved function with a U- 
shaped trend: at small masses the concentration first de- 
clines, reaches a minimum at log cr" 1 ~ 0.15, and then in- 
creases again at larger masses. The c — a relation is much 
narrower than the traditional c(M, 2) relation, but it is not 
exactly universal: there are some small dependences with 
redshift and cosmology. 

• The median concentration-<r(M2oo) relation can be ac- 
curately parameterized by eqs. (|14H16j) . This relation pro- 
vides an analytical model of the halo mass-concentration 
median relation c(M2oo) that reproduces all the relevant fea- 
tures, namely the decline of concentration with mass, and 
its flattening and upturn at high redshift /mass. 

• Our estimates for concentration of cluster-size ha- 
los are compatible with the recent observational results 
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and are substantially - a factor of 1.5 - larger than in 
iMaccio. Dutton. fc van den Bosch] l|2008h . 
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Figure 14. Effects of relaxation on halo concentration at rcdsliift 
z = 2. Top symbols show only relaxed halos. Concentrations of all 
halos are shown by bottom symbols. Relaxed halos are more con- 
centrated and, just as all halos, show an upturn at large masses. 
Open symbols show results from the Bolshoi simulation; filled 
symbols are for the Multidark simulation. Halos are binned by 
Vmax- The average halo mass A/200 for each velocity bin is shown 
on the horizontal axis. 



7 APPENDIX: EFFECTS OF HALO 
SELECTION 

Selection of halos has some impact on estimates of halo con- 
centrations. There are different ways to select halos and each 
selection condition has its own effect. One of the possibili- 
ties is to select quiet halos from a population of all halos. 
Halos can be selected by mass or by the maximum circular 
velocity V max , as we do in this paper. Halo radius may also 
be defined differently, which substantially changes the halo 
concentration. Here we consider and discuss the impact of 
different selection effects. 

The selection conditions used in the main part of this 
paper are: (1) all halos regardless of their relaxation status; 
(2) overdensity limit of 200 relative to the critical density 
in eq. ((2} defines the radius of a halo. In the main part of 
the paper we use notation c for concentrations defined in 
this way. In the Appendix we also use another definition 
of the virial radius. So, in order to distinguish these two 
definitions we will use C200 for the overdensity 200 definition 
and Cvir for the radius defined by the solution of the top-hat 
collapse model in t he AC DM cosmology as approximated by 
iBrvan Sz Normanl |T998). For the parameters of the Bolshoi 
and MultiDark simulations A v i r = 97 at z = 0. Relative to 
the matter density this corresponds to A = A/f2 m attcr = 
360. 

One of the interesting features in halo concentration 
is the upturn in concentrations at very large masses. Be- 
cause the largest halos also grow very fast, one may wonder 
whether the upturn is just a non-equilibrium feature. We 
address this issue by comparing concentrations of all ha- 
los with the concentrations of only relaxed halos at z — 2. 



Results presented in Figure [14] show that concentrations of 
relaxed halos are larger than for all halos by about 10% and 
more importantly that the upturn is also clearly present in 
relaxed halos. 

In Figure[l5]we investigate effects of selection by mass, 
effects of radius definition, and effects of relaxation at z = 0. 
On both panels the bottom curves are for all halos selected 
by mass and the top curves are for relaxed halos selected by 
Vmax. For the latter we show average mass for halos selected 
by velocities. Results for all halos selected by Vmax are in be- 
tween the two curves and are not shown t o avoid crowding . 
In agreement with previous results (e.g ., iNeto et"aL 1 120071 : 
M actio. Dutton, fc van den Boschl [20081 ), we find that con- 
centrations of relaxed halos are larger than those of all halos. 
For smaller masses M < IO 13 /i -1 M@ the two concentra- 
tions have nearly the same slopes and differ by ~ 7 — 10%. 
The difference increases at larger masses where the con- 
centration of relaxed halos practically stops declining at 
M w 10 15 /i _1 M Q . 

The following approximations provide fits for halo con- 
centrations at z = for masses M = 5 x I0 10 - 10 15 /i _1 M©: 
Relaxed halos selected by Vmax: 



C200 = 7.! 



/ M 2 , 



Vl0 12 ft -1 M e 



10.5 



Mvir 



1 + 0.2 



I + 0.15 



M 2l 



10 16 fe _1 M P 



Mvir 



1/2' 
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1/2-1 
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All halos selected by halo mass: 

M200 



C200 = 7.28 



10 12 /i _1 M P 



9.7 



M vir 



10 12 /i _1 M© 



(25) 
(26) 
(27) 



In summary, we find that different selection criteria 
work very differently: some are important and some are not: 

• Differences in concentrations due to selection of ha- 
los by mass or by Vmax are small at small masses: for 
M200 < IO 13 /i _1 M0 halos selected by V max have concentra- 
tions larger only by a factor 1.02—1.04. The differences some- 
what increase at larger masses: at M200 = 5 x 10 14 /i _1 Mq 
the difference is a factor of 1.08. 

• As expected, the overdensity threshold (in our case 
either 200p cr or the virial overdensity) has a large im- 
pact on halo concentration. We find that at z — the 
ratio of the concentrations is nearly independent on halo 
mass : c v ir / C 200 = 1.35. This is consistent with p revious 
results (e.g. iMaccio. Dutton. fc van den Bos ch 2008). Note 
that one may naively expect that there should be a mass- 
dependent correctlon_due the change in the overdensity 
threshold |Hu fc Kravtsovll2003h . However, this is only true 
for a toy model of an isolated halo with the NFW density 
profile. In reality, the situation is more complex and there 
are two effects: one due to the fact that the same halo is 
measured at different radii and another due to the differ- 
ence in the sets of halos: some of distinct halos defined by 
small R200 radius cease to be distinct and become subhalos 
for large virial radius. 
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Figure 15. Effects of halo selection. The left panel shows results for the overdensity 200 values; the right panel is for the virial overdensity. 
Bottom symbols and curves show concentrations for all halos selected by mass. The top symbols and curves are for relaxed halos selected 
by maximum circular velocity. Open symbols show results from the Bolshoi simulation; filled symbols are for the Multidark simulation. 
All halos selected by V max are in between the two sets of curves. 

• Halos with large masses are effected the most by a par- 
ticular choice of halo selection. For example, Figure [15] in- 
dicates that the concentration C200 of relaxed halos selected 
by V^ax is 20 percent larger than that of all halos selected 
by mass for M200 = 10 15 ft _1 M Q . 

This paper has been typeset from a TffjX/ I^TffjX file prepared 
by the author. 



